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^ ! 1. Introduction 

Short bursts with a duratio n of less than 2 s repre sent about 
7-H ; 25% of the BATSE sample dKouveliotou et alJI 1993I) but had to 
, wait until 200 5 (i.e. eight years after long bursts) to enter the 
afterglow era dGehrels et al.l 2005; Fo x et al.ll2005l iHiorth et al] 
2005; lBergeretal.ll2005l) . This is due to two reasons: (/) short 
bursts tend to e mit less photons becaus e of harder spectra and 
lower fiuences (Kouveliotou et al.ll 19931) . which makes their lo- 
calization more difficult; (ii) they have fainter afterglows, which 
are harder to detect. Following the discovery of the first after- 
glows, it appeared that the nature of the host galaxy, the loca- 
tion of the afterglow, and the absence of a supernova imprint 
in the visible light curve (even when the host is located at a 
redshift below 0.5) were indicativ e of progenitors tha t were dif- 
feren t from those of long bu rsts dGehrels et al.ll2005l iFox et al] 
120051: ISoderber g et al. 2006). Several short burts are clearly as- 
socia ted to elliptical galaxies (Blo om & Prochask a 2006; Bergerl 
2009) while others with accurate positions appear to have no 
coincident hosts, which clearly excludes progenitors belonging 
to the young pop ulation and favors merger scenarios involving 
compact objects {Naravan et al.lll992trM ochkov itch et al.I l993; 
iRuffert & Jankalll999UBelczvnski et al.ll2006l) ~ 

About 40% of the short bursts have no detectable afterglows 
after about 1000 s while the other 60% (ISakamoto & Gehrelsl 
2009) have long-lasting afterglows compar able to those of long 
bursts (see the review on short bursts by Nakar 2007 and ref- 
erences therein). If short bursts indeed result from the merging 
of two compact objects, the kick received when the black hole 



Context. GRB 080503, detected by Swift, belongs to the class of bursts whose prompt phase consists of an initial short spike followed 
by a longer soft tail. It did not show any transition to a regular afterglow at the end of the prompt emission but exhibited a surprising 
rebrightening after one day. 

Aims. We aim to explain this rebrightening with two different scenarios - refreshed shocks or a density clump in the circumburst 
medium - and two models for the origin of the afterglow, the standard one where it comes from the forward shock, and an alternative 
one where it results from a long-lived reverse shock. 

Methods. We computed afterglow light curves either using a single-zone approximation for the shocked region or a detailed multi- 
zone method that more accurately accounts for the compression of the material. 

Results. We find that in several of the considered cases the detailed model must be used to obtain a reliable description of the shock 
dynamics. The density clump scenario is not favored. We confirm previous results that the presence of the clump has little effect on 
the forward shock emission, except if the microphysics parameters evolve when the shock enters the clump. Moreover, we find that 
the rebrightening from the reverse shock is also too weak when it is calculated with the multi-zone method. On the other hand, in 
the refreshed-shock scenario both the forward and reverse shock models provide satisfactory fits of the data under some additional 
conditions on the distribution of the Lorentz factor in the ejecta and the beaming angle of the relativistic outflow. 

Key words. Gamma rays bursts: general; Gamma rays bursts: individual: GRB 080503; Shock waves; Radiation mechanisms: non 
thermal 



or neutron star components formed in past supernova explosions 
dHobbs et alJl2005l;lBelczvnski et alll2006l) can allow the system 
to reach the low-density outskirts of the host galaxy (or even to 
leave the galaxy) before coalescence occurs. This can naturally 
explains why some afterglows are so di m or have no coinc ident 
host (the observational data presented in lTroia et al.l l2008 show 
that the galactocentric offset of short bursts is on average much 
larger than for long bursts). 

The direct and simple connection between duration and 
progenitor class became fuzzier when it was found that in 
some bursts an initial short duration s pike is followed by a 
soft tail lasting several tens of seconds (|Barthelmv etail [2005; 

1120051: iNorris & Bonnellll2006l) . It was then sug- 
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gested (Zhan 1 120061: Gehr els et al.1 120061) to introduce a new 
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terminology that would distinguish type-I bursts resulting from 
mergers and type-II events coming from collapsars. In the ab- 
sence of a detected afterglow that can help to relate the burst to 
either the old or young stellar population, a vanishing spectral 
lag (for both the spike and the extended emission) has been pro- 
posed as an indicator for a type-I identification (Gehrel s~et al.l 
120061) . 

GRB 080503 belongs to the class of short bursts with ex- 
tended emission. The extended emission ended with a steep de- 
cay that was not immediately followed by a standard afterglow 
component. A peculiar feature in GRB 080503 is that after re- 
maining undetected for about one day, it showed a spectacu- 
lar rebrightening (both in X-rays and t he visible), which c ould 
be followed for five days in the visible. IPerlev et all (120091) de- 
scribed in great detail the multi-wavelength data they collected 
for this event and discussed different possibilities that could ac- 
count for the late rebrightening: (i) a delayed rise of the af- 
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terglow due to an extremely low density of the surrounding 
medium; (ii) the presence of a density clump in the burst en- 
vironment; (Hi) an off-axis jet that becomes visible when rel- 
ativistic beaming has been reduced by deceleration (see e.g. 
iGranot et al. 2002); (iv) a refreshed shock, when a slower part of 
the ejecta catch es up with the shock, a gain as a result of decel- 
eration (see e.g. lSari & M eszaros 2000); and finally (v) a "mini- 
supernova" from a small amount of ejec ted material powered by 
the decay of 56 Ni dLi & Paczvriskilll998l) . 

Case (0 imposes an external density below 10~ 6 cirT 3 , which 
seem s unreasonab ly low; case (Hi) implies a double-jet struc- 
ture (iGranotl 12005)" with one on-axis component producing the 
prompt emission but no visible afterglow (which can be possible 
only if the prompt phase has a very high efficiency) and the other 
one (off-axis) producing the delayed afterglow; case (v) can ac- 
count for the rebrightening in the visible, but not in X-rays. 

We therefore reconsider in this work the two most promising 
cases (zz) and (iv) in the context of the standar d model where the 
afterg l ow is produced by the forward shock dMeszaros & Reesl 
[19971 [SarietiD [l998) but also in the alternative one where 
it comes from a long-live d reverse shock dGenet et al.l l2007t 
lUhm & Belobo rodov 2007). The paper is organized as follows. 
We briefly summarize the observational data in Sect. |2]and list 
in Sect. [3] possible sources for the initial spike and extended 
emission. We constrain in Sect. |4]the energy released in these 
two components and discuss in Sect.[5]different ways to explain 
the rebrightening with a special emphasis on cases (ii) and (iv) 
above. Finally Sect.|6]is our conclusion. 

2. Short summary of the observational data 

2.1. Prompt emission 

The Swift-BAT light-curve of GRB 080503 consists of a short 
bright initial spike followed by a soft extended emission of re- 
sp ective durations f9o, sp ike = 0.32 + 0.0 7 s and f9o >ee = 170 + 40 
s (iMao et al.l l2008: Perl ev et al.1 12009). The fluence of the ex- 
tended emission from 5 to 140 s and between 15 and 150 keV 
was S,! e 5 ~ 150 = (1.86 + 0.14) 10~ 6 erg.cirT 2 while the fluence of 
the spike S was 30 times lower. The spectra of both the 
spike and the extended emission were fitted by single power- 
laws with respective photon indices 1.59 + 0.28 and 1.91 +0.12. 
The position of the initial spike in the duration-hardness diagram 
and the absence of any significant spectral lag (together with the 
absence of a candidate host galaxy directly at the burst loca- 
tion) make it consistent with a short (type I) burst classification, 
resulting from the merging of two compact objects. No spectral 
lag analysis could be performed on the extended emission, which 
was weaker and softer than the spike. 



2.2. Afterglow emission 

The afterglow of GRB 080503 was very peculiar. The prompt 
extended emission ended in X-rays with a steep decay phase of 
temporal index a = 2-4 (F(t) oc r a ), which is common to most 
long and short bursts. This decay did not show any transition 
to a "regular afterglow" and went below the detection limit in 
less than one hour. This behavior has been observed in about 
40% of the short burst population (ISakamoto & Gehrelsll2009l) 
but in GRB 080503 it covered nearly six orders of magnitude. 
In the visible, except for a single Gemini g band detection at 
0.05 day, the afterglow remained undetected until it exhibited 
a surprising late rebrightening (both in X-rays and the visible) 
starting at about one day after trigger. Following the peak of the 



rebrightening, the available optical data points (extending up to 
five days) and subseq uent upper limits show a steep decay of 
temporal index a ~ 2 dPerlev et al.ll2009l) . 

3. Origin of the different emission components 

The different temporal and spectral properties of the prompt ini- 
tial spike and extended emission indicate that they are produced 
by distinct parts of the outflow, possibly even with different dis- 
sipation or radiative mechanisms. The temporal structure of the 
extended emission, showing a short time-scale variability (with 
?var ~ 1 s), excludes the possibility of any conventional after- 
glow origin. Models of the central engine have been proposed, 
which are able to produce a relativistic outflow made of two 
distinct components with kinetic powers and temporal proper- 
ties similar to what is seen in short GRBs with extended emis- 
sion. For example, in compact binary progenitors, the extended 
emission coul d be caused by the fallback of m aterial, following 
coalescence (Rosswog 2007; Troi a et aLll2008l) . For a magnetar 
progenitor. iMetzger et al.l (120081) suggested that the initial spike 
is produced by accretion onto the protomagnetar from a small 
disk, while the extended emission comes from rotational energy 
extrac ted on a longer time scale. Finally, iBarkov & Pozanenkol 
(1201 ll) recently described a two-component jet model that could 
explain short GRBs both with and without extended emission, 
where a wide, short-lived jet is powered by vv annihilation and 
a narrow, long-lived one by the Blandford-Znajek mechanism. 

For the rest of this study it will be assumed that the outflow in 
GRB 080503 consisted of two main sub-components, responsi- 
ble for the initial spike and the extended emission, respectively, 
and that the afterglow emission is associated to the interaction 
of this structured outflow with the circumburst medium. The en- 
ergy content of each component can be estimated from the ob- 
served fluences (Sect. 0J. For the refreshed-shock scenario (see 
Sect. !5.2l below) their typical Lorentz factors are somewhat con- 
strained by the time of the rebrightening for a given value of the 
external density. 

4. Kinetic energy of the outflow 

To obtain the kinetic energy carried by the different parts of 
the outflow, one should start estimating the correction factor 
between the 15 - 150 keV and bolometric fluences C bo1 = 
S hol /S 15-150 for both components. Unfortunately, the shape of 
the spectrum is poorly constrained so that we will simply as- 
sume that 2 < C bo1 < 4. This is the range obtained with the 
simplifying assumption that the spectrum can be represented 
by a broken power-law of low and high-energy photon indices 
a = -1.5, f3 = -2.5, and peak energy between 20 and 300 keV 
(with C bo1 ~ 2 - 2.5 for E v between 20 and 100 keV and rising to 
4 at E p = 300 keV). From the fluence, we can express the total 
isotropic energy release in gamma-rays as a function of redshift 
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where £>l(z) is the luminosity distance. To finally obtain the ki- 
netic energy, one has to assume a radiative efficiency f m ^, defined 
as the fraction of the initial kinetic energy of the flow eventually 
converted to gamma-rays. The remaining energy at the end of 
the prompt phase is then given by 
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We adopted f m g » 0. 1 as a typical value It could be lower for in- 
ternal shocks dRees & Meszaroslll994t iDaigne & Mochkovitch| 
19981) or higher 



2001 



for ma g netic r econnection (Spruit et al. 
Drenkhahn & Spruitl 120021: iGiannios & Spruitl 120061 



McKinnev & Uzdensky 2010) or modified photospheric emis- 
sion dRees & Meszarosl 120051; iBel oborodov 120101) . We did not 



consider scenarii where the radiative efficiency would be very 
different for the spike and extended emission even if this possi- 
bility cannot be excluded a priori. Because the redshift of GRB 
080503 is not known, we adopted z = 0.5 as a "typical" value for 
a type-I burst for the different examples considered in Sect. [5] 
This yields fi K ,i S0 - C h s f x 1.1 10 52 erg and C*^ x 3 1 50 erg 
for the extended emission and spike. The dominant uncertainties 
on these energies clearly come from the unknown radiative effi- 
ciency and distance of the burst. We briefly discuss below how 
our results are affected when assuming a different redshift or a 
different radiative efficiency. 



that is injected in the population of accelerated relativistic elec- 
trons (power-law distribution with a slope —p) and in the am- 

l ii 

plified magnetic field - follow the prescription e e = e B , which 
results from the acceleration process of el ectrons moving t oward 
current filaments in the shocked material (lMedvedevl2006l) . This 
assumption simplifies the discussion but is not critical for the 
general conclusions of our study. 

We did not try to fit the initial steep decay in X-rays because 
it is generally interpreted as the high-latitude emission ending 
the prompt phase and not as a true afterglow component. In that 
respect, it is not clear if the optical data point at ~ 0.05 day 
should be associated to the high-latitude emission or already be- 
longs to the afterglow. We assumed that it is of afterglow origin 
(the most constraining option) and imposed that the simulated 
light curve goes through it. This leads to some specific conse- 
quences, mainly for the reverse shock model (see discussion in 
Sect. ESS. 



5. Modeling the afterglow of GRB 080503 

5. 1 . Forward and long-lived reverse shocks 

We considered two different mechanisms that can explain GRB 
afterglows. The first one corresponds to the standard picture 
where the afterglow results from the forward shock propagat- 
ing in the external medium, following the i nitial energy depo- 
sition by the cen tral engine (Sari et al.|[l998b. The second one 
was p roposed by iGenet et al.l (I2007I) and lUhm & Be loborodov 
(2007) to account for some of the unexpected features revealed 
by Swift observations of the early afterglow. It considers that the 
forward shock is still present but radiatively inefficient and that 
the emission comes from the reverse shock that sweeps back into 
the ejecta as it is decelerated. The reverse shock is long-lived be- 
cause it is supposed that the ejecta contains a tail of material with 
low Lorentz factor (possibly going down to T = 1). 

We performed the afterglow simulations using two different 
methods to model the shocked mat erial. In the first o ne it is rep- 
resented by one single zone as in ISari et aL | (119981) : the phys- 
ical conditions just behind the shock are applied to the whole 
shocked material. At any given time, all shocked electrons are 
considered as a single population, injected at the shock with 
a power-law energy distribution. Then the corresponding syn- 
chrotron spectrum can be calculated, taking into account the ef- 
fect of electron cooling over a dynamical timescale. The sec- 
ond method is more accurate, consideri ng separately the e volu- 
tion of each elementary shocked shell (Beloborodov 2005J) ex- 
cept for the pressure, which is uniform throughout the whole 
shocked ejecta. The electron population (power-law distribution) 
and magnetic field of each newly shocked shell are computed 
taking into account the corresponding shock physical conditions 
and microphysics parameters. Then each electron population is 
followed individually during the whole evolution, starting from 
the moment of injection, and taking into account radiative and 
adiabatic cooling. The evolution of the magnetic field - assum- 
ing that the toroidal component is dominant - is estimated using 
the flux conservation condition. Furthermore, it was checked that 
the magnetic energy density never exceeds equipartition. 

Finally, we made a few more assumptions to somewhat re- 
strict the parameter space of the study. We adopted a uniform ex- 
ternal medium of low density because GRB 080503 was proba- 
bly a type-I burst, resulting from the coalescence of two compact 
objects in a binary system at the periphery of its host galaxy. We 
also assumed that the redistribution microphysics parameters e e 
and es - respectively the fraction of the shock dissipated energy 



5.2. Refreshed shocks 

One way to explain the late rebrightening is to consider that the 
forward or reverse shocks have been refreshed by a late supply of 
energy dRees & Meszaroslll998HSari & Meszarosll2000l) . This is 
possible if the initial short duration spike in the burst profile was 
produced by a "fast" relativistic outflow (of Lorentz factor r sp jk e ) 
while the extended emission came from "slower" material with 
Tee < Tspike- Then, at early times, only the fast part of the flow 
is decelerated and contributes to the afterglow. When the slower 
part is finally able to catch up, energy is added to the shocks and 
the emission is rebrightened. 

5.2.1. Forward shock model 

In the standard forward shock model the lack of any detectable 
afterglow component before one day imposes severe constraints 
on either the density of the external medium or the values of the 
microphysics parameters. Fixing e e and e# to the com monly used 
value s 0.1 and 0.01 implies to take n < 10~ 6 cm" 3 dPerlev et al.l 
2009). This very low density would likely correspond to the in- 
tergalactic medium, which might be consistent with the absence 
of a candidate host galaxy down to a visual magnitude of 28.5. 
We preferred to adopt a less extreme value n = 10~ 3 crrr 3 , more 
typical of the interstellar m edium at the outskirts of a galaxy 
(see e.g. Steidel et al. 2010). Then, decreasing the microphysics 

1 /2 

parameters to e e = e B = 0.05 becomes necessary to remain 
consistent with the data. 

To obtain a rebrightning at one day we adopted F ee = 20 
and F sp ike = 300. The outflow lasts for a total duration of 100 s 
(Is for the the spike and 99 s for the tail). We injected a kinetic 
energy E^n = 7 10 50 erg in the spike and 50 times more in the 
tail. It can be seen that the results, shown in Fig.[TJ are consistent 
with the available data and upper limits except possibly after the 
peak of the rebrightening where the decline of the synthetic light 
curve is not steep enough. This can be corrected if a jet break oc- 
curs close to the peak, which is possible if the jet opening angle 
6>jet is on the order of 1 /r ee 0.05 rad. This beaming angle is 
somewhat smaller than the values usually inferred from obser- 
vations of short b urst afterglows (see e.g. iBurrows et al.l 12006: 
iGrupe et all 120061) or suggested by simulations of com pact bi- 
nary mergers (see e.g.|Rosswog & Ramirez-Ruizll2002l) . 

An example of a light curve with a jet break is shown 
in Fig. [TJ assuming that the jet has an opening angle of 3.4° 
(0.06 rad) and is seen on-axis. A detailed study of the jet- 
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break properties is beyond the scope of this paper and we 
therefore did not consider the case of an off-axis observer 
and neglected the lateral spreading of the jet, expected to 

become important when r < \J 6\e-t. Detailed hydrodynami- 

cal studies (see e.g. Granot 20071; Zhang & MacFadven| [2009; 



Ivan Eerten & M acFadvenf 201 U Lvutikovf 2011b tend to show, 
however, that as long as the outflow remains relativistic, the jet- 
break is more caused by the "missing" sideways emitting mate- 
rial than by jet angular spreading. 

We finally checked how our results are affected if different 
model parameters are adopted. If the density n of the external 
medium is increased or decreased, similar light curves can be 
obtained by changing the Lorentz factors (to still achieve the 
rebrightning at one day) and the microphysics parameters (to 
recover the observed flux). For example, increasing the density 
to n = 0.1 crrT 3 requires r ee =* 10 (keeping F^e = 300) and 
e e = e B 2 = 0.02. Conversely, with n = 10~ 5 cm 4 , F ee — 35 and 

1 /2 

e e — e B = 0.08 are needed. 

If the kinetic energy of the outflow is increased (resp. de- 
creased) because the radiative efficiency / ra( j is lower (resp. 
higher) or the redshift higher (resp. lower), light curves agree- 
ing with the data can again be obtained by increasing (resp. de- 
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creasing) F ee and decreasing (resp. increasing )e e = e B '\ Also 
note that the spread of the Lorentz factor 5r ee around F ee at 
the end of the prompt phase has to be limited to ensure that 
the slower material is able to catch up in a sufficiently short 
time to produce an effective rebrightening. In the case shown 
in Fig. [1] we have 5r ee /T ee = but we have checked from 
the numerical simulation that acceptable solutions can be ob- 
tained as long as £r e e/r e e $ 0.2. This configuration is for ex- 
ample naturally expected after an internal shock phase where 
fast and slow parts of the flow collide, resulting in a shocked re- 
gion with a nearly uniform L orentz factor distribution (see e.g. 
iDaigne & Mochkovitchl2000l) . 

5.2.2. Long-lived reverse shock model 

If the afterglow is produced by the reverse shock, similar good 
fits of the data can be obtained. Fig.|2]shows an example of syn- 



thetic light curves for E\ 
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spike 
kin 



= 7 1 50 erg and Eg, = 30 E 



spike 
J kin ' 



= 0.16, p - 2.5 in the shocked ejecta and n = 10 
cm -3 . The microphysics parameters have to be higher than in 
the forward shock case because the reverse shock is dynami- 
cally less efficient, which requires a higher radiative efficiency 
to obtain the same observed fluxes. The Lorentz factor distribu- 
tion in the ejecta is also slightly different to guarantee that the 
light curve (i) goes through the optical point at 0.05 day and (ii) 
decays steeply after the peak. No jet break has to be invoked 
here because the decay rate (in contrast to what happens in the 
forward shock model) depends on the distribution of energy as a 
function of the Lorentz factor in the ejecta. 

Again, if the external density and kinetic energy of the flow 
are varied, satisfactory fits of the data can be recovered by 
slightly adjusting the Lorentz factor and microphysics param- 
eters. 



5.3. Density clump in the external medium 

We now investigate the possibility that the rebrightening is 
caused by the encounter of the decelerating ejecta with a den- 
sity clump in the external medium. For illustration, we adopted 
a simple distribution of the Lorentz factor that linearly decreases 



with injection time from 300 to 2 so that, in the absence of the 
density clump, afterglow light curves from either the forward or 
reverse shocks would be smooth and regular (we have checked 
that the exact shape of the low Lorentz factor tail is not crucial in 
this scenario). To model the clump, we assumed that the circum- 
burst medium is uniform (with n = 10~ 3 cm 4 ) up to 1.7 10 18 
cm (0.55 pc) and that the density then rises linearly to n — 1 
cm -3 over a distance of 10 18 cm (0.32 pc). The ejecta is strongly 
decelerated after entering the high-density region and we find 
that the forward shock is still inside the clump at the end of the 
calculation (at t \, s = 8 days). 

5.3.1. Forward shock model 



As iNakar & Granotl (I2007I) showed by coupling their hydrody- 
namical calculation to a detailed radiative code, a density clump 
in the external medium has little effect on the forward shock 
emission. Therefore a clump cannot produce the rebrightening 
in GRB 080503. In the simple case where the shocked medium 
is represented by a single zone, the effect of the clump is barely 
visible. With the detailed multi-zone model a stronger rebright- 
ening is found, because the effects of the compression resulting 
from the deceleration of the flow are better described, but even 
in this case the calculated flux remains nearly one order of mag- 
nitude below the data. 

Fig. [3] illustrates these results and confirms that the forward 
shock emission does not strongly react to the density clump. 
Even if, from an hydrodynamical point of view, the forward 
shock is sensitive to the clump, the observed synchrotron emis- 
sion is only moderately affected because the increase in up- 
stream density is nearly counterbalanced by the decrease of 
Lorentz factor in the shocked material. Of course, spectral ef- 
fects complicat e the picture, but the e ssence of the result remains 
the same (see lNakar & Granot]|2007l for details). 

5.3.2. Possible evolution of the microphysics parameters 

In view of the many uncertainties in the physics of collisionless 
shocks it is often assumed for simplicity, as we did so far, that 
the microphysics redistribution parameters e e and stay con- 
stant during the whole afterglow evolution. However, particle- 
in- cell simulations of accelera tion in collisionless shocks (see 
e.g. lSironi & Spitkovsky|201ll) do not show any evidence of uni- 
versal values of the parameters. If e e or/and eg are allowed to 
change during afterglow evolution, the problem of the forward 
shock encountering a density clump can be reconsidered, now 
with the possibility of a sudden increase of radiative efficiency 
triggered by the jump in external density. 

Fig. |4] shows the resulting light curves when the micro- 
physics parameters e e and eg of the forward shock are increased 

1 /2 

at the density clump. If the prescription e e = e B is maintained, 
no satisfactory solution can be found in the simple model where 
the shocked medium is represented by one single zone. In this 
case, the optical frequency lies between the injection and cool- 
ing frequencies (v, < v opt < v c ) while the X-ray frequency 
satisfies vx > v c so that the visible and X-ray flux densities 
depend on the microphysic s parameters in the following way 
dPanaitescu & K umar 200^) 



/y,o P t « ef e B 4 and / v>x oc e'J L e B 4 



P- 1 . 



(3) 



With the prescription e e = e B we obtain / v opt oc e e 2 and / v> x oc 

3p-4 

e e 2 The optical flux is therefore much more sensitive than the 
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Fig. 1. Refreshed shocks: forward shock model. Left panel: initial distribution of the Lorentz factor (lower part) and kinetic power 
(upper part) in the flow as a function of inject ion time t m \ . Right p anel: synthetic light curves at 2 eV (black, dotted line) and 10 keV 
(gray, dotted line) compared to the data from lPerlev et al.l (120091) . The kinetic energies injected in the spike and extended emission 
components are E^ e = 7 IO 50 erg and E£? = 50 E^ e . We adopt e e = e^ 1 = 5 10~ 2 , p = 2.5 in the shocked external medium 
together with n = 10~ 3 cirT 1 and z = 0.5. The steeper thin lines at late times correspond to a conical jet (seen on-axis) of opening 
angle 6j et = 0.06 rad. 
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Fig. 4. Density clump: forward shock model with varying mi- 
crophysics parameters. Synthetic light-curves at 2 eV (black) 
and 10 keV (grey) when the microphysics parameters of the for- 
ward shock are changed at the density clump: e e is increased by 
a factor 5 from an initial value of 10~ 2 , keeping the prescription 

1 II 

e e — e B . The Lorentz factor, injected kinetic power in the out- 
flow, and the density distribution in the external medium are the 
same as in Fig. [3] The dotted and dashed lines correspond to the 
single and multi-zone models for the shocked region. The same 
simulations with the addition of a jet-break (assuming 9j et = 0.08 
rd) are also shown. 



X-ray flux to a change of the microphysics parameters and a 
simultaneous fit of the data in both energy bands is not possible. 
A simple solution to this problem is to change e e alone, keeping 
€b constant. In this case, increasing e e by a factor of 25 (from 
0.01 to 0.25) is required to reproduce the rebrightening in both 
the X-ray and visible ranges. 

In the more detailed model with a multi-zone shocked re- 
gion, the situation is different. In the shells that contribute most 
to the emission, we find that both v opt and vx are larger than v c 
and therefore / vopt and / v x depend in the same way on the mi- 
crophysics parameters. It is then possible to achieve a satisfac- 
tory solution (dashed lines in Fig. |4j that keeps the prescription 
e e = e l J 2 (with e e increased by a factor 5). As in Sect. I5.2.1l we 
introduce a jet break (now assuming 6>j e t = 0.08 rd) to account for 
the decay of the optical flux following the peak of the rebright- 
ening. Notice that the decay is steeper here (compare Fig. [T]and 
Fig. HJ owing to the rapid decrease of the Lorentz factor inside 
the clump. 

5.3.3. Long-lived reverse shock model 

With the simple one-zone model, the reverse shock emission is 
found to be much more sensitive to the density clump than the 
forward shock emission. Indeed, when the ejecta starts to be de- 
celerated, its bulk Lorentz factor suddenly decreases and slow 
shells from the tail material pile up at a high temporal rate and 
with a strong contrast in Lorentz factor. These two combined 
effects lead to a sharp rise of the flux from the reverse shock. 
Synthetic light curves showing a satisfactory agreement with the 
data are shown in Fig. [3] 

However, the detailed multi-zone model gives different re- 
sults, where the rebrightening is dimmer and cannot fit the data. 
The main reason is that the higher contrast in Lorentz factor, 
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Fig. 2. Refreshed shocks: long-lived reverse shock model. Left panel: initial distribution of the Lorentz factor (lower part) and 
kinetic power (upper part) in the flow as a function of injection time f; n j . Right panel: synthetic light-curves at 2 eV (black, dotted 
line) and 10 keV (grey, dotted line) together with the data. The kinetic energies in the spike and extended emission components are 
£y * e = 7 1 50 erg and E^r = 30 E^ e . The density of the external medium, redshift and slope p of the electron distribution are the 
same as in Fig.Q] The adopted microphysics parameters are e e = e l J 2 = 0.16 in the shocked ejecta. 
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Fig. 3. Density clump: forward and long-lived reverse shock models. Left panel: initial distribution of the Lorentz factor (lower 
part) and kinetic power (upper part) in the flow as a function of injection time t m y Middle and right panels: synthetic light-curves at 
2 eV (black) and 10 keV (grey) for the forward and reverse shocks, using either the simple one-zone (dotted lines) or the detailed 
multi-zone (dashed lines) model. The kinetic energies in the spike and extended emission components are E^ e = 7 1 50 erg and 

^kin = 30 E^ e . The adopted microphysics parameters are e e - e l J 2 = 10~ 2 (forward shock, middle panel) and e e = e^ 1 - 0.07 
(reverse shock, right panel). See text for the prescription adopted for the density clump. 



leading to a higher specific dissipated energy, now only concerns 
the freshly shocked shells, while in the single zone model it is ap- 
plied to the whole shocked region. This example (as well as the 
one already discussed in Sect. I5.3.2I I shows that using a detailed 
description of the shocked material can be crucial when dealing 
with complex scenarios (i.e. not the standard picture where the 
blast-wave propagates in a smooth external medium, with con- 
stant microphysics parameters )Q. 



It is still possible to fit the data by increasing the mi- 
crophysics parameters during the propagation in the clump. 
However, this seems less natural than for the forward shock 
(Sect. 15.3.2b because the upstream density of the reverse shock 
does not change. On the other hand, a modification of the mi- 
crophysics parameters could still be due to the sudden increase 
in the reverse shock Lorentz factor triggered by the clump en- 
counter. 



1 In the refreshed-shock scenario (Sect. 15.21 where the dynamics is 
simpler, the single and multi-zone models give comparable results. 
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6. Conclusion 

GRB 080503 belongs to the special group of short bursts where 
an initial bright spike is followed by an extended soft emission of 
much longer duration. It did not show a transition to a standard 
afterglow after the steep decay observed in X-rays at the end of 
the extended emission. This behavior has been observed previ- 
ously in short bursts, but GRB 080503 was peculiar because it 
exhibited a spectacular rebrightening after one day, both in X- 
rays and the visible. The presence of the extended emission pre- 
vents one from classifying GRB 080503 on the basis of duration 
only, but the lack of any candidate host galaxy at the location of 
the burst and the vanishing spectral lag of the spike component 
are consistent with its identification as a type-I event resulting 
from the coalescence of a binary system consisting of two com- 
pact objects. 

From its formation to the coalescence, the system can mi- 
grate to the external regions of the host galaxy allowing the burst 
to occur in a very low density environment, accounting for the 
initial lack of a detectable afterglow. To explain the late rebright- 
ening, we considered two possible scenarios - refreshed shocks 
from a late supply of energy or a density clump in the circum- 
burst medium - and two models for the origin of the afterglow, 
the standard one where it comes from the forward shock and the 
alternative one where it is made by a long-lived reverse shock. 

In the refreshed-shock scenario we supposed that the ini- 
tial spike was produced by fast moving material (we adopted 
T = 300) while the one making the soft tail was slower (T ~ 20). 
Initially, only the spike material is decelerated and contributes 
to the afterglow until the tail material is eventually able to catch 
up, which produces the rebrightening. Both the forward and re- 
verse shock models provide satisfactory fits of the data under the 
condition that the material making the tail has a limited spread 
in Lorentz factor 6T/T < 0.2. This allows the rise time of the 
rebrightening to be sufficiently short. This condition might be 
satisfied from the beginning but can also result from a previ- 
ous sequence of internal shocks that has smoothed most of the 
fluctuations of the Lorentz factor initially present in the flow. 
In addition, a jet break is required in the forward shock case to 
reproduce the steep decline that follows the rebrightening . This 
implies that the jet should be beamed within an opening angle of 
3 - 5°, which appears somewhat smaller than the values usually 
preferred for type-I bursts. In the long-lived reverse shock model 
a jet break is not necessary because the shape of the light curve 
now depends on the energy distribution in the ejecta, which can 
be adjusted to fit the data. 

In the scenario where a density clump is present in the burst 
environment, the rebrightening resulting from th e forward shock 
is wea k, in agreement with the previous work of Naka r & GranotI 
(2007). We performed the calculation in two ways: first with 
a simple method where the shocked material was represented 
by one single zone, then using a more detailed, multi-zone ap- 
proach. The impact of the clump was barely visible in the first 
case. The rebrightening was larger in the second one but still re- 
mained nearly one order of magnitude below the data. We then 
considered the possibility that the shock microphysics might 

change inside the clump. We found that by increasing e e by a 

l ii 

factor of five (and with the prescription that e e = e B ) it was 
possible to fit the data with the multi-zone model under the ad- 
ditional condition to have a jet break at about 2-3 days (corre- 
sponding to a jet opening angle < 5°). With the simplified model 
the results were more extreme, imposing to increase e e alone by 
a very large factor of 25. 



If the afterglow is made by the reverse shock, the effect of 
the clump is strong with the simple model. It is however much 
reduced with the detailed model and the observed rebrighten- 
ing cannot be reproduced, the synthetic light curve lying nearly 
one order of magnitude below the observed one. It appears that 
only the multi-zone approach provides a proper description of 
the compression resulting from the encounter with the density 
barrier. Conversely, in the refreshed-shock scenario the simple 
and detailed models give comparable results. 

From the different possibilities we considered, which could 
explain the late rebrightening in GRB 080503, several appear 
compatible with the data, but none is clearly favored. The 
refreshed-shock scenario may seem more natural because the 
initial spike and extended emission probably correspond to dif- 
ferent phases of central engine activity. It is not unreasonable to 
suppose that the material responsible for the extended emission 
had a lower Lorentz factor, as required by the refreshed-shock 
scenario. Then, both the forward and reverse shock models lead 
to satisfactory fits of the X-ray and visible light curves, if two 
conditions on the Lorentz factor distribution and jet opening an- 
gle (see above) are satisfied. The density clump scenario does 
not seem able to account for the rebrightening if the afterglow 
is made by the reverse shock. The conclusion is the same with 
the forward shock, except if the microphysics parameters are al- 
lowed to change when the shock enters the clump. 
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